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A quadratic potential in a light cone QCD inspired model 
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Abstract. The general equation from previous work is specialized to a quadratic potential V(r) = —a + 
\fr 2 acting in the space of spherically symmetric S wave functions. The fine and hyperfine interaction 
creates then a position dependent mass m(r) in the effective kinetic energy of the associated Schrodinger 
equation. The results are compared with the available experimental and theoretical spectral data on the 
7r and p. Solving the eigenvalue problem within the usual oscillator approach induces a certain amount of 
arbitrariness. Despite of this, the agreement with experimental data is within the experimental error and 
better than other calculations, including Godfrey and Isgur |^| and Baldicchi and Prosperi |1()|. The short 
coming can be removed easily in more elaborate work. 

PACS. ll.10.Ef - 12.38.Aw - 12.38.Lg - 12.39.-x 



1 The S-state Hamiltonian 

For spherically symmetric S states the previously derived 
Hamiltonian reduces in Fourier approximation to [T||2"||!J| 
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There are no more interactions than the central potential, 
the hyperfine, the kinetic, and the Darwin interaction, but 
also no less. For s-states the total spin squared is a good 
quantum number S 2 — [(ci + 02)/2] 2 = 5(5 + 1), thus 



o-i cr 2 = 25(5+1) -3 
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1, triplet, 
0, singlet. 
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Shaping notation, the Hamiltonian is written as 
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The non locality of the Hamiltonian resides in the position 
dependent mass 
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To solve this Hamiltonian, one must go on a computer. 

The Hamiltonian in Eq.JSJ looks like a conventional 
instant form Hamiltonian as obtained by quantizing the 
system at equal usual time. But it must be emphasized 
that it continues to be a genuine front form or light cone 
Hamiltonian 5 , derived from the latter by a series of exact 
unitary transformations 



Because it is shorter, cti(T2 is kept explicit in the equations 
as an abbreviation for Eq.J2J. With a quadratic potential, 
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the spring constant is /, the Hamiltonian Q becomes the 
non-local Schrodinger equation 
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2 The model Hamiltonian and its parameters 

In this first round, I try to avoid to go on the computer 
as far as possible, by the following reason. According to 
renormalization theory, the renormalization group invari- 
ants (parameters) must be determined from experiment. 
This is a strongly non linear problem. In order to get a first 
and rough estimate, the Hamiltonian is simplified here un- 
til it has a form which is amenable to analytical solution. 
Therefore, all in-tractable terms in the above will be re- 
placed here by mean values and related to the experimen- 
tally accessible mean square radius (r 2 ) [H]. In effect, the 
substitution 
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Table 1. Model parameters in GeV. Note: [/*] = 1. The heavy 



Table 2. Model parameters in GeV. Note: [/* 



flavor masses are determined by the sing 


det triplet difference. 
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0.4768 | 0.5899 | 1.7646 | 5.1555 


0.6291 | 3.4035 || 


|| 4 || 0.2465 | 0.2648 | 0.2657 | 0.3481 


|| 0.5412 | 1.7598 || 









is the only true assumption in the present model. I con- 
sider thus the model Hamiltonian, 
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Its eigenvalues are 

E n = — a + oj £o + oj rj n + caia 2 , 



(8) 

(9) 
(10) 

(11) 
(12) 



with £o 



| and i] n — In. The invariant mass squares 



M 2 = (?77l + 777 2 ) 

+ 2(mi + m 2 )( - a + £qlu + rj n uj + caia 2 ) , (13) 

are then related to experiment. 

For equal masses mi = m 2 = m, the model has the 3 
parameters m, / and a. One thus needs 3 empirical data 
to determine them. I choose: 
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The spectrum is labeled self explanatory by the flavor 
composition M n = Mdu,tn or M n = M du ,sn, for singlets or 
triplets, respectively. The triple chosen in Eq. l|14|l exposes 
a certain asymmetry. The excited p is chosen since its ex- 
perimental limit of error is very much smaller than the 
one for the corresponding tt state. Only its ground state 
mass is known very accurately, i.e. m w + = 139.57018 ± 
0.00035 MeV. In the present work only the first 4 digits 
are used. For equal masses, the above abbreviations be- 
come 



c 



= 1 



/ 

2m 2 

a + 3c , 

7(/(r 2 )/2-q) 
8m 



(15) 



The experiment defines 2 certainly positive differences: 

(16) 
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A third one can be constructed by the observation that 
^X 2 - |y 2 - Mj fi s0 = 4ma - 4m 2 . Keeping in mind 

that oj 2 = 2//m, one can remove trivial kinematic factors 
and define 3 experimental quantities B, C and D by 
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Substituting / = mC 2 and ma = B 2 
7 
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a quadratic equation with the solution 

D 2 + IB 2 
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Having m, the / and a are calculated from l|17|) . 

The position dependent mass changes the relation be- 
tween the mean square radius and it experimental value. 
Therefore, I introduce a fudge factor /* according to 



(r 2 ) = (fl 2 (r% 



(19) 



Since all mesons have about the same size by order of 
magnitude, this number is kept universal. The fudge fac- 
tor is introduced here to account, in some global fashion, 
for the tremendous simplification introduced by replacing 
Eq.© with (JHJ. Some large scale variations are compiled 
in Table |3 The mass spectra including the ground states 
vary very little with the fudge factor. Any variations would 
show up the fastest for the high excitations. For this rea- 
son, the masses for 77 = 4 are included in the table. I 
do not understand this insensitivity from a mathematical 
or numerical point of view. The major effect of /* is the 
ease by which one can change the quark mass. A value of 
/* ~ 40 leads to the 20 MeV for the quark mass quoted 
in ^U]. Here /* = 4 is chosen. The other components of 
the spectrum, both for the p and the tt, are then obtained 
for free. They will be compiled in Tabled] below. 



Table 3. Dependence on the fudge factor. 
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0.850 


0.367 


0.935 


6.065 


2.1648 


2.2929 


2.0 


0.477 


1.165 


0.629 


3.403 


2.1648 


2.2929 


4.0 


0.247 


4.358 


0.541 


1.760 


2.1648 


2.2929 


8.0 


0.124 


17.131 


0.709 


0.888 


2.1648 


2.2929 


16.0 


0.062 


68.220 


1.228 


0.445 


2.1648 


2.2929 


32.0 


0.031 


272.58 


2.362 


0.223 


2.1648 


2.2929 
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Table 4. S wave spectra in GeV for light unflavored mesons. 



Table 5. S wave spectra in GeV for strange mesons. 



n 


% Singlets 


7T+(ud) 


n 


3 Si Triplets 


p+(ud) 




Experiment 


Theory 




Experiment 1 


Theory 


1 


0.1396(0) 


0.1396 


1 


0.7685(6) 


0.7685 






0.1396 2 






0.7711 2 






0.150 3 






0.769 3 






0.497 4 






0.846 4 
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1.300(100) 


1.2550 


2 


1.465(25) 


1.4650 






1.2650 2 






1.4650 2 






1.300 3 






0.769 3 






1.326 4 






1.461 4 
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1.795(10) 


1.7694 


3 


1.700(20)" 


1.9240 






1.7950 2 






1.9230 2 






1.880 3 






2.000 3 






1.815 4 






1.916 4 


4 




2.1648 
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2.150(17) 


2.2929 
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2.2912 2 



Hagiwara etal [7], Zhou and Pauli |8]. 
3 Godfrey and Isgur 4 Baldicchi and Prosperi |1U| (a) 
"Could be a D state [TT|. 



In principle, one could determine the heavier quark 
masses analytically from the hyperfme splittings: 
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The so obtained results are however, not very reasonable, 
see Tabled The experimental numbers are insufficiently 
accurate. Therefore, I determine them numerically from 
the singlets M uStS o, M„g jS o and M u ^ s o and compile them 
in Table |2 



3 Results and Discussion 

Unflavored light mesons. The results for the ir-p sys- 
tem are compiled in Table 0] The experimental points are 
taken from from Hagiwara et al 0. It is no surprise that 
theory and experiment coincide for the 7r + , the p + and 
the p + (1450), because these data have been used to de- 
termine the parameters. The remarkable thing is that one 
can perform such a fit at all, The model reproduces the 
huge mass of the excited pion within the error limit. This 
solves the long standing puzzle, why a physical system can 
have a first excited state with a ten times larger mass. 

The remaining three calculated masses of the ir-p sec- 
tor agree with experiment almost within the error bars. 
The model underestimates the second 7r-excitation by only 
26 MeV. The second excitation of the p is overestimated 
by a comparatively large 224 MeV, but the experiment for 
the p + (3 3 Si) needs confirmation. The third excitation of 
the p + (4 3 Si) is overestimated by 224 MeV. 

The table includes also a comparison with other the- 
oretical calculations. It includes the results from a recent 
oscillator model [Hj- Their model is even simpler than the 



n Experiment 1 Theory n Experiment 1 Theory 



'So Singlets K+(us) 



0.493677(16) 



1.460" 



1.830 a 



0.4937 

0.6048 2 

0.47 3 

1.3732 

1.5480 2 

1.45 3 

1.8782 

2.1040 2 

2.02 3 

2.2736 



'Si Triplets K* + (us) 



0.89166(26) 



1.629(27) i ' 



0.0.8718 

0.8917 2 

0.90 3 

1.5498 

1.6808 2 

1.58 3 

2.0109 

2.6242 2 

2. II 3 

2.3845 



1 Hagiwara etal 7 , 2 Zhou and Pauli 8 . 

3 Godfrey and Isgur 0, 

"To be confirmed; b J p not confirmed. 



present one: it works with a hyperfme splitting, only, but 
suppresses the mechanism of a position dependent mass. 
Despite this, the results of [8] coincide practically with the 
present ones. I have included also the results from the pi- 
oneering work of Godfrey and Isgur [3] as a prototype of a 
phenomenological model, and from a recent advanced cal- 
culation by Baldicchi and Prosperi (TJJ- Neither of these 
models have much in common with the present one. They 
fail to reproduce the pion, this mystery particle of QCD. 

Strange mesons. The S wave K + and K* + spectra 
are given in Table The mass of the singlet ground state 
is used to determine the mass parameter m s . Except the 
ground states, the experiments carry many ambiguities 
about the quantum number assignment for K and K* 
mesons. The model prediction for the triplet ground state 
underestimates the experimental value by 20 MeV. Both 
the first and the second excited state of K (2 1 So and 3 1 So) 
are not confirmed. Another unconfirmed resonance with 
mass 1.629 ±0.027 GeV lying between 2 1 So and 3 1 So was 
assigned to be a singlet K. Apparently there is no position 
for it in the K spectrum if it is an S wave state. However, 
according to its mass, it might well be the first excited 
state of K* (2 1 So). Taken the numbers in the table, the 
discrepancies are 88 and 69 MeV for the singlet and triplet 
n = 2 states, respectively. The second excited state of the 
K (2 1 S ) differs by only 48 MeV, but the datum needs 
confirmation. 

Heavy mesons. The S wave uc, ub, sc, sb and cb 
meson spectra are given in Table No excitations were 
observed for these mesons. 

The uc singlet is used to determine the mass param- 
eter m c . Its ground state mass (5°) therefore coincides 
with experiment. The model overestimates the mass of 
the triplet D*° (1%) by about 50 MeV. — The ub singlet 
is used to determine the mass parameter mt. Its ground 
state mass (B + ) therefore coincides with experiment. The 
model overestimates the mass of the triplet B* + (l 1 So) by 
16 MeV only. 

No data in the sc mesons are used to determine model 
parameters. Model and experiment differ by 27 and 40 
MeV for singlet and triplet, respectively. 
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Table 6. Ground state masses in GeV for heavy mesons. Table 7. The predicted S spectrum in GeV for heavy mesons. 



n Experiment 1 Theory n Experiment 1 Theory 



% Singlet D"(uc) 



1.8645(5) 



1.9961 
1.9224 2 

i qq3 



% Singlet B+(ub) 



5.2790(5) 



5.2790 

5.2965 2 

5.31 3 



x So Singlet D~(sc) 



1.9685(6) 



1.9961 
2.0201 2 



% Singlet B»(sb) 



5.3696(24) 



5.3976 

5.3739 2 

5.35 3 



'So Singlet B+(cb) 



6.4(4) 



6.5077 
6.27 3 



5 Si Triplet D* u (uc) 



2.0067(5) 



2.0718 

2.0067 2 

2.04 3 



% Triplet B* + (ub) 



5.3250(6) 



5.3085 

5.3250 2 

5.37 3 



3 Si Triplet D* a ~ (sc) 



2.1124(7) 



2.0718 

2.0655 2 

2.13 3 



'Si Triplet B* u (sb) 



5.4166(35) 



5.4214 

5.3885 2 

5.45 3 



3 Si Triplet B*+(cb) 



6.5157 

6.3458 2 

6.34 3 



Hagiwara etal 7 , 2 Zhou and Pauli 8 . 
3 Godfrey and Isgur 
a To be confirmed; b J p not confirmed. 



n L S 


D u 


(uc) 


n^Si 


zr u 


(uc) 


1 

2 
3 
4 


1.8645 
2.5013 
3.0061 
3.4375 


1 

2 
3 
4 


1.9594 
2.5728 
3.0658 
3.4899 


r/So 


B+ 


(ub) 


n 3 Si 


B*+ 


(ub) 


1 

2 
3 
4 


5.2790 
5.8473 
6.3651 
6.8438 


1 

2 
3 
4 


5.3085 
5.8739 
6.3896 
6.8666 


r/So 


D7 


(sc) 


n^Si 


D* s - 


(sc) 


1 

2 
3 
4 


1.9961 
2.5837 
3.0605 
3.4724 


1 

2 
3 
4 


2.0718 
2.6426 
3.1104 
3.5165 


n So 


B» 


(sb) 


n Si 


BT 


(sb) 


1 

2 
3 
4 


5.3976 
5.9163 
6.3930 
6.8365 


1 

2 
3 
4 


5.4214 
5.9380 
6.4131 
6.8553 


n'So 


Bt 


(cb) 


n a Si 


B* c + 


(cb) 


1 

2 
3 
4 


6.5077 
6.8447 
7.1659 
7.4733 


1 

2 
3 
4 


6.5157 
6.8523 
7.1731 
7.4802 



Model and experiment differ by 27 and 4 MeV for sin- 
glet and triplet, respectively. 

Model and experiment agree for the singlet. The triplet 
data are unknown. 

The model prediction for the complete spectrum are 
compiled in Table [7| for easy reference. 

The flavor diagonal mesons like ss, sc or bb may not 
be calculated in the model, see [H]- 

4 Conclusions 

The agreement between the present simple model with an 
oscillator potential and the experiment is generally good. 
These are good news, since harmonic interactions are easy 
to work with in many body problems. The present ap- 
proach will be useful for considering baryons and nuclei. 

With the 4 mass parameters of the up/down, strange, 
charm and bottom quarks, the model has only 2 addi- 
tional 2 parameters for the harmonic oscillator potential. 
In principle, the fudge factor must be counted as parame- 
ter, but as seen above, the choice of the up/down and the 
fudge factor is strongly coupled. 

The 6 canonical parameters of the model generate a 
reasonably good agreement with the 21 data points avail- 
able. 

Note that renormalized gauge field theory has also 
4+1+1 parameters: The 4 flavor quark masses, the strong 
coupling constant a s , and the renormalization scale A. Of 
course, they can be mapped into each other PEE]- 

Once one has determined the parameters in such a first 
guess, one should relax the model assumption, Eq.JJJ), and 



work with the full non local model, with a position depen- 
dent mass. For this one has to go back to the computer 
and perform the necessary fine tunings of the parameters. 
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